Polariton parametric photoluminescence in spatially inhomogeneous systems 
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A general theory of polariton parametric photoluminescence in spatially inhomogeneous systems 
is developed. The quantum Langevin equations are solved in a generalized Bogoliubov de Gennes 
approximation. We apply the formalism to the specific case of a disordered microcavity. In this 
case, we numerically solve the equations for the coherent emission and the photoluminescence. We 
describe the effect of the exciton and photon disorder on the photoluminescence pattern exhibited 
in momentum space, finding a good agreement with the experimental observations. 
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Polariton photoluminescence is the paradigm for gen- 
eration of quantum correlated states of electronic excita- 
tions in semiconductor devices. A pair of quantum corre- 
lated polaritons are generated via the scattering process 
2kp ks -\- ki^ i.e. two polaritons in the pumped mode 
kp are scattered into one polariton in the "signal" mode 
ks and one polariton in the "idler" mode ki [1]. Since 
all the k-modes satisfying the energy-momentum conser- 
vation rules are admitted as final modes, the resulting 
two-polariton state is an entangled state. Experimen- 
tal evidence of the quantum nature of the correlations 
between the signal and the idler polaritons has been re- 
ported [1, 0] . Polaritons are also advocated as a vector of 
quantum correlations, in particular allowing the genera- 
tion and the control of long-range quantum correlations 
that are stored in other systems, for example providing 
the ideal system for realizing spin coupling within a very 
long range [4]. This possible application is very promis- 
ing, because of the recent engineering of new microstruc- 
tures able to trap polaritons in spatially inhomogeneous 
regions @, Q . These polariton traps are the ideal can- 
didates for designing polariton quantum devices. This 
scenario calls for a general theory of parametric photolu- 
minescence in presence of spatial inhomogeneity. In ad- 
dition, the irnportant role of structural disorder has been 
highlighted 0, H, @], in particular in relation to para- 
metric photoluminescence 0, Isl • 1^ been shown that 
the photoluminescence from idler modes is suppressed 
because it is very sensitive to disorder on the exciton 
component. 

In this letter, we develop a general formalism for de- 
scribing parametric photoluminescence in presence of ar- 
bitrary spatial inhomogeneity. The theory is based on 
the Bogoliubov de Gennes approach, typically used for 
describing Bose-Einstein condensation (BEG) in inhomo- 
geneous systems. Bogoliubov theory assumes the pres- 
ence of a large coherent field, that is treated classically, 
resulting in a linearized Hamiltonian for the quantum 
fluctuation field. In the present case, the classical Bo- 
goliubov field is induced by the pump, while the residual 
quantum fluctuation field describes the actual paramet- 
ric photoluminescence. Our theory is formulated in the 



basis of the exciton and photon fields and includes the 
linear exciton-photon coupling and the typical exciton 
non-linearities [10]. After deriving the general theory, 
we present the predictions for a disordered system and 
discuss the comparison with experiments ^ . 

We treat the exciton and the photon fields as 
Bose fields and we consider the effective Hamiltonian 



H = Hq -\- Hr -\- Hx -\- Hs 



(1) 



where = J dr^Ur)[-{n^ /2mx)\/^ + /7,(r)]^,(r) + 
/ (ir^J(r)[ec(V) + /7c(r)]^c(r) is the unperturbed 
term, Hr = Qr J (ir[^J.(r)^c(r) + h.c] is the 
term describing the exciton-photon coupling, Hx = 
l^a; /(ir^J.(r)^J.(r)^^(r)^2;(r) is the effective 2-body 
exciton interaction term, modeling both Goulomb inter- 
action and the effect of Pauli exclusion on electrons and 
holes 0, and Hs = Vs J dr[^{r)^{r)^x{r)^x{r) + 
h.c] is the term modeling the saturation of the exciton 
oscillator strength [10]. Here we neglect the momentum 
dependence of the scattering matrix elements Vx and Vs^ 
because for our present purposes we will consider only 
scattering processes between states with small momen- 
tum. 

We consider the system evolving under a continuous 
monochromatic optical pump, F{r,t) = e~*'^p^F^(r), and 
we assume that the exciton and photon fields decay in 
time with the rates jx and 7c, respectively. In this non- 
equilibrium regime, the two fields are written as |13| 



^x(c)(r,t) 



^x(c)(r) + V^a:(c)(r,t) 



(2) 



i.e. as the sum of a classical term {'^x{c){'^^t)) = 
e~*^p^$a,(c)(r), describing the coherent field generated by 
the pump and evolving accordingly with the pump fre- 
quency cjp, and a Bose field tljx(c){'^^t)^ describing the 
fluctuations. By averaging the Heisenberg equations of 
motion of x{c){'^^t)^ we obtain the two coupled equa- 
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tions 



2mr 



+2t;,Re{$;(r)$e(r)})*x(r) 
{nn + Vs\Mr) |2) $,(r), 



(3) 



hojp<i>c{r) = ec(V) +/7c(r) -i7c $c(r) 



+ {nn + Vs\<i>Ar) |2)$.(r)+F0(r), (4) 

defining tiie spatial siiape ^x,c{'^) of tiie coiierent fields. 

By adopting an input-output formalism [3] and by 
linearizing the time evolution of the excitation field, we 
derive the quantum Langevin equation for the fiuctua- 
tions 11 5l 



(5) 



where ^ = {ipx^ipl^ipc^'^l)^ is the four-component 
fiuctuation field, while f = (/x, /J, /c^ /c) defines the 
fiuctuation field, with correlations (/^(r, t))f = and 
{fdr,t)fj{r',t'))t = 27r^^5(,^5{r - r')S{t - t') [15]. The 
matrix M has the Bogoliubov de Gennes form [16] 



M = 



-(S-f)* 

\ -(Sf2^)* 



Sxx 
12 

-Tx - ilx 

S^xc 







^12 
-^R 


-fc - i-ic ) 



where 



+ U^ijc) - fkOp + 2v^ I $x(r) 



2mx 

+ 4«,Re{$:(r)$e(r)}, 
S?f(r) = t;,$2(r)+2t;,$,(r)$e(r), 
nnir) = nR + 2vs I $x(r) I' , 
Efl(r) = v^'^lir). 



(6) 



(7) 
(8) 
(9) 
(10) 



We notice that the interaction terms in Eqs. (|3|), (|4|) and 
(j6j) are written within the Bogoliubov limit, i.e. they only 
depend on the coherent fields, because we have assumed 
that the incoherent populations are vanishingly small. 
This is always true in the regime of parametric photolu- 
minescence that we are considering in the present work. 
However, in the same spirit as for symmetry breaking 
theories of BEG [17], our symmetry breaking approach 
allows to systematically include the higher order interac- 
tion processes involving excitations, and thus to treat a 
regime of higher density [l|, [l^ . 

By Fourier transforming Eq. ([5|) in the frequency do- 
main 



hcj^ir, co) = M^(r, co) + f (r, co) , (11) 
cy resolved exciton (photon) 



we have access to the frequency resolved exciton (photon) 
spatial density na,(c)(i*7^) ~ 



using the correlation properties of the fiuctuation field 
and after some algebra, we can write the densities in the 
useful form 

n,(,)(r,c^) = ^ / dr^|(r\/|MjV,l(3))f r,, (12) 

where T = (0, 27r7^, 0, 27r7c)^ and = M - Tiujl is 
evaluated in the basis |(r, j)), where |r) spans the 2D real 
space, while the label j = 1, ...,4 refers to the bloc form 
of M, Eq. (|6]). In this work, we focus on the regime of 
spontaneous photoluminescence and we neglect the den- 
sity of excitations in the expression of the correlations 
^^ 

The excitation fields (k, uj) expressed in the in-plane 
momentum space obey an equation formally equivalent 
to Eq. (pTj) and the densities in momentum space are 
given by 



(k,a;) 



4 

E 

1=1 



dk' 



{k',?|Mji|k,l(3)) 



(13) 



Full information about the photoluminescence emission 
is obtained from Eqs. (p!2|) and (fT3|) . that can be nu- 
merically solved for any shape of the external potentials 
Ux,c{'^) and of the exciting pump F^(r). The form of Eqs. 
(p!2]) and (fT3|) is particularly advantageous for numerical 
computations, because of the efficiency of the algorithms 
for the treatment of sparse matrices. 

We now apply the model to the specific situation in 
which the external potentials Ux,c{'^) describe the struc- 
tural disorder naturally present in the polariton system 
[9]. We adopt parameters modeling a typical GaAs 
microcavity [uj. In particular, we assume radiative 
linewidths j^^^ = 0.06 meV and j'^^^ = 0.13 meV, and 
we assume that the external pump have frequency hcOp = 
— 1.2 meV with respect to the heavy- hole exciton ground 
state and in-plane momentum kp = (1.73,0) /im~^. The 
photon disorder potential is assumed to be Gauss corre- 
lated in space ^ 



(C^c(r)C/c(r')) 



(14) 



with correlation amplitude <Jc = 0.1 meV and correlation 
length = 7 /im. These values result into an inhomoge- 
neous broadening AEc ~ 50 /ieV for the cavity photon 
resonance. The exciton disorder is assumed with a short 
correlation (white noise) on our computation grid and 
its amplitude is fixed in order to give an inhomogeneous 
broadening AEx ~ 0.75 meV for the exciton resonance 
[9]. 

The solution of Eqs. ^ and ^ is plotted in Fig. [H 
In Fig. [TJa) we show the square modulus of the resulting 
coherent photon field |$c(i*)P in real space. The photon 
disorder is responsible for the localization of the field over 
a distance of tens of microns, while the exciton disorder 
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FIG. 1: (Color online) (a) Square modulus of the coherent 
field l^clr)!^ in the real space, in linear color scale, (b) 
Square modulus of the coherent field |$c(k)|^ in the in-plane 
momentum space, in logarithmic color scale. Notice the spot 
at (1.73,0) /im~^ corresponding to the incident pump field 
(intensity off scale). 



of the eight-shape is due to the finite linewidths), clearly 
appears. We call "Signal" emission that for kx < kp^, 
while "idler" emission that for kx > kp^. In the dis- 
ordered case, Fig. [2] (a), the momentum selection rule 
is lifted and the pattern is broadened and fragmented. 
At large momenta, the effect is dramatic, resulting into 
a complete spreading out of the idler emission over the 
excitonic states. The vanishing of the idler resonance is 
due to the strong effect of the exciton disorder in the 
large-momenta region of the polariton dispersion, where 
the very large effective mass makes the states sensitive 
to the disorder-induced broadening of the exciton res- 
onance. This result agrees with the experimental ob- 
servations, because a clear evidence of the eight-pattern 
has been reported only for a sample specially designed 



to minimize exciton disorder As expected from the 
linear response theory, at small momenta the effect of ex- 
citon disorder is less evident 0. However, in Fig. [2](a) 



is responsible for the additional short-length fluctuations 
of the amplitude of the field. The square modulus of 
the coherent photon field |$c(k)p in the {kx,ky)-spdice 
is shown in Fig. [TJ^b). Here, the typical elastic ring of 
Rayleigh scattered polaritons is clearly visible p. 9, 20l, 

El. 
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FIG. 2: Pattern in the momentum space of the photolumi- 
nescence for the disordered microcavity (a) and for a uniform 
microcavity (b), plotted in logarithmic color map. 

We then solve Eq. ([13]) and compute, in the momen- 
tum space, the parametric photoluminescence intensity, 
which corresponds to the frequency-integrated photon 
density /pi,(k) = J dconcCk^oo). In Fig. [2] we compare 
the resulting photoluminescence intensity (in panel (a)) 
with the one obtained for the uniform system (in panel 
(b)), i.e. in the absence of photon and exciton disorder. 
In the uniform case, a eight-shaped photoluminescence 
pattern, due to the momentum selection rules (the width 
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FIG. 3: Polariton population np{kx,ky) in presence of disor- 
der (a) and for a uniform system (b). The color map is in 
linear grey scale. 

we also notice the presence of speckles, both in the idler 
and in the signal emission. The qualitative difference be- 
tween their distribution and the one of the speckles in 
the linear response (Fig. [1]) is probably related to the 
non-linear nature of the parametric process. 

Although speckles are hardly visible in photolumines- 
cence, they are an indication of the dramatic effect of the 
exciton disorder on the polariton population np(k, cj) = 
nc(k, cj) + naj(k, cj), i.e. the quantity which is more rele- 
vant for future applications in the generation of quantum 
correlated pairs of polaritons. This quantity is considered 
in Fig. [3l where we compare the resulting frequency- 
integrated polariton population (panel (a)) with the one 
obtained for the uniform system (panel (b)). The ef- 
fect of disorder is even more striking. The spots in the 
idler region indicate the fragmentation in the momentum 
space of the population distribution. This fragmentation 
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also reflects on the population distribution at small mo- 
menta, which is broadened with respect to the homoge- 
neous case. Notice that in the uniform case, because of 
the signal-idler symmetry, the relative amplitude of the 
frequency integrated populations of the signal and idler 
modes depends on the ratio between linewidths, via the 
simple relation 7kn(k) = 72kp-k^(2kp — k). In the disor- 
dered case, because of the localization of the modes, this 
relation is no longer valid. 
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FIG. 4: Frequency integrated polariton population np(|k — 
kp|) along the cut shown in the inset. We compare the results 
obtained for a uniform system (dashed line), in presence of 
photon disorder only (dot-dashed) and in presence of both 
exciton and photon disorder (solid). 

Fragmentation is even more evident in Fig. [H where we 
show the frequency integrated polariton population along 
one specific cut in the momentum plane (the dashed line 
shown in the inset of the figure). In the absence of disor- 
der, the distribution is narrow around the signal and idler 
modes and the population is larger in the idler than in 
the signal, because of the longer radiative lifetime. When 
photon disorder is introduced, a broadening appears, cor- 
responding to the softening of the momentum selection 
rule, but this effect is relatively small and the signal-idler 
pair is clearly visible. In the presence of exciton dis- 
order, the population in the idler region is significantly 
fragmented in momentum space. Correspondingly, the 
signal resonance is broadened. Indeed, the parametric 
process mixes the high energy exciton-like states and the 
low- moment a ones. Notice that the same effect can be 
seen in the idler when only photon disorder is present 
(dot-dashed line). We have checked that, while the dis- 
appearance of the eight-shape of the photoluminescence 
pattern and the broadening of the population distribu- 
tion in the signal can be qualitatively reproduced by in- 
troducing inhomogeneous linewidths 0], the fragmenta- 
tion of the polariton population, due to the presence of 
localized states, can be predicted only by a full treatment 
of the spatial disorder. 

In conclusion, we have developed a general formal- 
ism for describing the spontaneous parametric photolu- 



minescence in a spatially inhomogeneous polariton sys- 
tem. The range of application of this formalism is very 
large, allowing the description of disordered systems and 
of the recently engineered polariton traps [l, 0] . In the 
present work we have applied the formalism to study the 
appearance of the eight-shaped pattern in the photolu- 
minescence of a disordered microcavity. We have shown 
that the disordered-induced fragmentation of the popu- 
lation distribution in the momentum space is responsible 
for a strong suppression and broadening of both the idler 
and the signal resonances. The important effect on the 
signal population is a consequence of its parametric cou- 
pling with the idler. Our formalism is the tool required 
for designing structures for the generation of quantum 
correlated polariton states. 

We are grateful to I. Carusotto for enlightening discus- 
sions. 
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